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Abstract. Exotic Ni isotopes have been measured at the R 3 B-LAND setup at GSI in 
Darmstadt, using Coulomb excitation in inverse kinematics at beam energies around 500 MeV/u. 
As the experimental setup allows kinematically complete measurements, the excitation energy 
was reconstructed using the invariant mass method. The GDR and additional low-lying strength 
have been observed in 68 Ni, the latter exhausting 4.1(1.9)% of the El energy-weighted sum rule. 
Also, the branching ratio for the non-statistical decay of the excited 68 Ni nuclei was measured 
and amounts to 24(4)%. 
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1. Introduction 

Unstable nuclei have been shown to exhibit several exotic features that cannot be observed 
in their stable counterparts. In neutron-rich nuclei, one of these features is the appearance 
of electric-dipole strength at energies near the neutron separation threshold. Located below 
the well-known Giant Dipole Resonance (GDR), which has already been studied extensively 
in stable species [I], this new low-lying El strength distribution is often denoted as Pygmy 
Dipole Resonance (PDR), referring to its much smaller photoabsorption strength compared to 
that of the GDR. Despite its small El strength, the PDR provides an experimental access to 
the equation-of-state (EOS) of asymmetric nuclear matter, and thus a link to the neutron skin 
evolution [2J, as well as about the abundance distribution of heavy elements in the universe. 
Indeed, these two phenomena arise from the behavior of the EOS of nuclear matter, more 
precisely of its symmetry-energy term. By constraining the various parameters of the EOS 
using measurements of the PDR, not only the global description of nuclear matter under various 
conditions (asymmetry and density) can be improved, but also the size, structure and chemical 
composition of neutron stars can be evaluated in a more realistic manner. The description of 
the abundance of heavy nuclei, on the other hand, can be directly influenced by the spectral 
shape of the low-lying El strength [3], affecting the particle (neutron or proton) capture and 
photodissociation reactions in various nucleosynthesis scenarios. 

One of the first PDR observations was realized at the LAND setup by studying the Coulomb 
excitation (j*,xn) of neutron-rich Sn isotopes [3]. The PDR was measured in 130 Sn and 132 Sn 
at energies around 10 MeV and exhausted 7(3)% and 4(3)% of the El energy-weighted sum 
rule, respectively. The comparison of these results with RPA calculations [5] served to constrain 
the symmetry energy and the symmetry-energy pressure, as well as the neutron-skin thicknesses 
of the studied Sn isotopes and of 208 Pb. More recently, the low-lying El strength in 68 Ni was 
measured with (7*, 7') virtual photon scattering j6], revealing 5% of the El energy- weighted 
sum-rule strength under the assumption of a branching ratio for the direct 7-decay of ~ 4% in 
the PDR region. From these results, Carbone et al. [7] derived symmetry-energy parameters 
that are in good agreement with those previously obtained from 130,132 Sn. 

The scarce experimental data on low-lying El strength in exotic nuclei leaves many important 
questions related to the PDR open. For instance, is there a systematic dependence of the PDR 
energy, width or strength on variables such as the mass, charge or isospin asymmetry? Also, 
does the decay of the PDR significantly differ from that of the GDR, giving, e.g., enhanced 
direct 7 or direct neutron decay branching ratios? What is the degree of collectivity of the 
PDR, and is it's character rather isovector or isoscalar? All these questions require systematic 
studies of the low-lying El strength to be carried out over various mass and isospin-asymmetry 
ranges, in which quantities, such as the differential photoabsorption cross section or the dipole 
polarizability, are measured and can then be compared to theory. 

We present here an experiment in which the El strength has been measured in neutron-rich 
Ni isotopes by means of electromagnetic excitation in inverse and complete kinematics. After the 
description of the R 3 B-LAND experimental setup and of the measurement principle, preliminary 
results for 68 Ni will be presented. 

2. Experimental setup 

The R 3 B-LAND setup is located at GSI Helmholtzzentrum fur Schwerionenforschung GmbH in 
Darmstadt, Germany, and allows the investigation of reactions in inverse kinematics involving 
very short-lived nuclei. Stable primary ion beams are prepared and accelerated to relativistic 
energies in the UNILAC linear accelerator combined with the SIS- 18 heavy- ion synchrotron. 
The production of unstable secondary beams takes place in the fragment separator FRS [8], 
consisting of a production target and of a series of dipole magnets allowing for the in-flight 
separation of the produced species. In the present case of the investigation of 68 Ni, a stable 86 Kr 



Figure 1. (Color online) R 3 B-LAND setup for the Ni experiment at GSI, Darmstadt. See text 
for description. 

primary beam was accelerated to approximately 650 MeV/u and directed onto a 4.2 g/cm 2 Be 
target, generating multiple secondary species by projectile fragmentation at the entrance of the 
FRS. Within two distinct settings of the FRS, the Ni isotopic chain was covered from A=67 to 
A=72, with fragment energies of the order of 500 MeV/u. The exotic secondary beams contained 
a mixture of several nuclei of interest delivered simultaneously, therefore requiring an eventwise 
particle identification, making use of three measurements at the FRS: a precise measurement 
of the Bp value for each particle, its velocity from a time-of-flight measurement and its charge, 
obtained from the energy loss. 

The R 3 B-LAND setup was arranged in the present case to allow a kinematically complete 
measurement of electromagnetic excitation reactions. Due to the high beam velocities of the 
order of /3 ~ 0.75, all the reaction products are strongly forward-focused. In most cases, this 
effect is beneficial, since almost full solid-angle coverage can be achieved with reasonably-sized 
detectors. In other cases, however, the kinematical boost induces an additional difficulty for the 
detection and data analysis, such as for the Doppler shift in the photon decay of excited reaction 
fragments. Figure [T] shows the R 3 B-LAND setup used for the present experiment on neutron-rich 
Ni isotopes. The secondary beam enters the setup from the lower right side of the figure, and 
passes through a series of tracking and identification detectors. Two silicon pin-diodes (PSP1/2) 
provide the energy-loss and position measurements, while the plastic scintillator detector POS 
generates a signal used as minimum-bias trigger and as timing reference for all other detectors. 
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Figure 2. (Color online) Virtual photon spectra for 68 Ni at 500 MeV/u on a na *Pb target. 
Virtual photon numbers for El (solid red line), E2 (dotted green line) and Ml (dashed blue 
line) photons. 



The reaction target is located directly in front of the Csl gamma detector, a 2tt barrel-type 
detector consisting of 144 crystals arranged in 12 azimuthal and 12 polar segments, detecting 
photons emitted by the excited beam fragments. The fine segmentation was designed to provide 
a relatively small Doppler broadening, thus increasing the energy resolution of the detector. 
Heavy charged particles, both non-reacted incoming species as well as reaction fragments, 
are strongly forward- focused and are detected by a third pin-diode (PSP3), providing a first 
energy- loss measurement of the reaction products and a position measurement. These charged 
particles are then deflected in the ALADIN dipole magnet into a fragment branch equipped with 
three plastic scintillating fiber detectors (GFI1-3) [9] and a plastic scintillator time-of-flight wall 
(TFW). These detectors provide position, time-of-flight and energy- loss measurements, from 
which the fragment mass can be reconstructed. Neutrons evaporated during the decay process 
of the excited nuclei also have very small laboratory angles, allowing them to reach the Large 
Area Neutron Detector LAND [10], located at 0° with respect to the incoming beam axis and 
approximately 15 m downstream of the target. 

Coulomb excitation is performed preferentially using a high-Z target, such as na *Pb in the 
present experiment. Figure [2] displays the number of virtual photons for as a function of their 
energy [11], representing the availability of a given excitation energy. This particular case was 
calculated for a 68 Ni beam impinging on a na *Pb target at 500 MeV/u, with a minimum impact 
parameter of 12.98 fm, obtained with the parameterization of Benesh, Cook and Vary [12] based 
on the nuclear radii of the collision partners. The nature of the Coulomb excitation process 
prevents us from knowing the excitation energy a priori. Therefore, it can only be reconstructed 
using the invariant-mass technique, in which the four-momenta of all species participating in 
the reaction must be measured. The reconstructed excitation energy can be expressed as: 



E * = /H m i + Hlj m i m j (! - Pifij cos + E lySum - m proj (1) 

V * 

where the indices i and j both run over all outgoing species and E 1)Sum is the total energy carried 
by photons. It should be noted that equation [1] is only valid in this form when all quantities are 



expressed in the incoming projectile's rest frame. While choosing an appropriate target nucleus 
to maximize the probability of electromagnetic excitation, a non-negligible probability for nuclear 
excitation remains, as well as a reaction-like signature induced by background interactions 
of the beam particles. As the differential cross section measured with the Pb target can be 
approximated as a linear combination of Coulomb, nuclear and background contributions, each 
nucleus of interest must be measured on three different targets, namely nat Pb, nat Q and an empty 
target. These three measurements allow a separation of the various contributions, providing 
clean Coulomb excitation differential cross sections. 

3. Evaluation of experimental response 

The comparison of theoretical differential cross sections with the experimental distributions 
requires an accurate understanding of the response function induced by the detectors of the 
setup and by the algorithms used to treat the measured data. Mainly two detectors of the setup 
are responsible for the response, namely LAND and the Csl detector. The neutron reconstruction 
algorithm of LAND processes the observed hits of each event into neutron hits using a series of 
conditions on the time and position variables. Due to this algorithm, there is a non-negligible 
probability of not detecting the proper number of neutrons. Additionally, the energy-dependent 
efficiency and acceptance must also be taken into account at this level. In the case of the 
Csl detector, the main response contribution lies in the limited calorimetric properties of the 
detector, thus shifting the real photon energies to much lower values. All of the effects cited 
above are evaluated quantitatively using a proper decay simulation for the nucleus of interest 
linked to a response simulation, taking into account all detector-related effects. While the 
Csl response is mainly simulated using GEANT4 [13] within the R3BROOT framework |14j . 
LAND relies on a different simulation principle. Instead of simulating completely the response 
of the neutron events with a simulation package {e.g., GEANT4), the neutrons from the event 
generator are projected onto the current LAND configuration, taking its distance to the target 
and eventual inactive sub- modules into account. For each neutron hitting LAND, a previously 
measured neutron hit distribution, including all secondary particles, is retrieved from a database 
and shifted to the simulated interaction point. All neutron hits for a given event are processed 
in this manner, before the complete hit distribution is evaluated by the neutron reconstruction 
algorithm. The neutron hit distributions in the database have been previously measured with 
LAND and consist of real 1-neutron events from a deuteron break-up experiment. This procedure 
allows us to be independent of any simulation package for the description of the interaction of 
high-energy neutrons with matter. 

Based on this simulation of the detector response, we are able to compare any theoretical 
input with measured data. This requires the convolution of the theoretical input according to 
the following expression for, e.g., the excitation energy: 




where E g is the generated excitation energy, E r is the reconstructed excitation energy, i and j 
are the number of generated and detected neutrons, respectively, 1Z is the response matrix, B 
is the branching ratio matrix and M is the mass separation matrix. The observable E of Eq. [2] 
can be exchanged freely for any other observable of interest. While the matrices 1Z and B are 
obtained from the statistical decay simulation, the mass separation matrix M is defined by the 
experimental conditions. This matrix takes into account the eventual limited mass separation 
capability of the setup, therefore mixing data from reaction channels involving neighboring 
fragment masses with the channel of interest. In the present case, the mass resolution obtained in 
the experiment is sufficient to exclude contaminants by mass cuts, such that the mass separation 
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Figure 3. (Color online) Response matrices for the 68 Ni(7*,2n) 66 Ni channel, relating the 
reconstructed excitation energy E r (left) and reconstructed neutron kinetic energy Efa nT (right) 
with the generated excitation energy E g . 
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The matrix B (i, E g ) consists of the branching ratio for channel A ~X.("y* ,ii\) A ~ l Y as a function of 
the generated excitation energy E g (or any other observable used in Eq. [2]) . 

The description of the response matrix lZ(i,j,E g ,E r ) is somewhat more complicated: it 
relates two observables, one generated and one reconstructed, for a given generated channel 
" 4 X(7*,m)" 4_l Y and for a given number of reconstructed neutrons j. Since in the present case the 
mass separation matrix does not contain off-diagonal elements, we only need response matrices 
where i = j, since all other combinations of i and j are rejected. Figure [3] shows two response 
matrices for the 68 Ni(7*,2n) 66 Ni channel (i = j = 2) for the reconstructed excitation energy 
E r (left) and reconstructed neutron kinetic energy Ekin.r (right) as a function of the generated 
excitation energy E g . The response matrix for the excitation energy nicely illustrates the strong 
distortion induced by the Csl detector. An ideal setup, consisting of fully efficient detectors 
and without any distortions due to the response, one would expect a strong linear correlation 
lying on the diagonal E r = E g line of the response matrix. This is indeed the case for small 
excitation energies (from S2n. up to ~15.5 MeV). At larger energies, however, the correlation 
on the diagonal line disappears in favor of a second correlation band with an offset to lower E r 
values. At even higher energies, more such correlation bands appear. This is due to the limited 
calorimetric efficiency of the Csl detector, as soon as photon emission is involved in the decay of 
the excited nucleus. In the case of 7-cascades (e.g., above 22 MeV in the left part of Fig. [3]), at 
least two secondary bands are visible, each one representing combinations of missed individual 
photons. 

Due to the complicated structure of the response matrices, a deconvolution of the measured 
data with the inverted response is not feasible. In order to compare a theoretical input with the 
experimental data, an iterative scheme based on the convolution of the theoretical input must 
be used. At the beginning of each iteration, a trial input, usually representing the shape of the 
expected El strength, is convoluted according to Eq. [21 producing the reconstructed differential 



cross section based on the input distribution. In the second step, the convoluted trial function 
is compared to the experimental data by calculating the x 2 value. The last step consists of 
an adjustment of the trial input in order to minimize the \ 2 value. The iterative procedure 
continues until a x 2 minimum is reached, at which point one considers that the trial input 
provides a possible deconvoluted description of the experimental data. 



4. Neutron decay of 68 Ni 

In the present experiment, the measured El strength was limited to excitation energies above 
the neutron threshold of the nucleus of interest. Therefore, several decay channels are open, 
such as the direct photon, direct particle and compound nucleus decays. When studying giant 
resonances in heavy nuclei, it is usually safe to assume that the compound nucleus decay is 
the dominant channel. Statistical models are capable of describing quantitatively the kinetic- 
energy distributions of particles emitted by an evaporation-like decay of the compound nucleus. 
In the case of one-neutron emission from a nucleus A, these distributions follow Maxwell 
distributions [T5l: 
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where a c is the neutron capture cross section on nucleus A-l and T is the nuclear temperature 
of the daughter nucleus. Since Eq. [3] is only valid for decays involving one neutron, Le Couteur 
and Lang [16] provide an expression valid for n > 1 emitted neutrons: 
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where I m yy and T* m j^T. The substantial difference between expressions 0] and [5] resides 
in the pre-factor of the exponential term. One therefore expects a pre-factor ~ Ekin in the In 
decay channel, and ~ ^ Ekin for all other neutron decay channels. 

The quantitative understanding of the experimental neutron kinetic-energy distributions 
requires a simulation of the decay of the compound nucleus coupled to a simulation of the 
detector response, as described in section The generation of neutron energies in this process 
is, however, difficult, since for a given excitation energy, the number of emitted neutrons is 
often not known a priori. Figure U] shows graphically how the algorithm chooses the Maxwell 
pre-factor. Only if the first candidate for the neutron evaporation is being processed and only 
if the excitation energy of the parent nucleus is located between the In and 2n thresholds the 
Ekin pre-factor is chosen. In all other cases, yJEki n is used. In order to calculate the kinetic- 
energy distributions, additional information concerning T is required. The temperature of a 
nucleus is usually defined as T = y/U/a, where U is the internal excitation energy and a is the 
level-density parameter of the nucleus. Evaluating a is not straightforward, however, several 
parameterizations exist. In the present case, we chose to use ABRABLA |17[ll8j to evaluate the 
nuclear temperature profiles for various nuclei of interest. Figure [5] presents the temperature of 
66 ~ 68 Ni as a function of excitation energy calculated in this manner. The curves are directly used 
by the statistical decay algorithm, allowing it to directly generate neutron energies according to 
Eq. land ED 

One can now attempt to describe quantitatively the energy-differential cross section of the 
evaporated neutrons in the 68 Ni(7*,n) 67 Ni channel, as shown in Fig. [6l The black points 
represent the measured data, obtained from the time-of-flight measurement between the target 
and LAND, as well as the position of the reconstructed neutron in the latter detector. As 
described in section [2J the differential Coulomb-excitation cross section is obtained by measuring 
the reaction on a nat Pb, a nat C and an empty target in order to disentangle the electromagnetic, 
nuclear and background components. While the background subtraction is straightforward and is 
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Figure 4. (Color online) Visualiza- 
tion of the algorithm used to gener- 
ate neutron kinetic energies. The pre- 
factor of the Maxwellian distribution is 
shown for the n th neutron being pro- 
cessed by the event generator. 
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Figure 5. (Color online) Temperature 
profiles for 68 Ni (solid red line), 67 Ni (dashed 
green line) and 66 Ni (dotted blue line) nuclei. 
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Figure 6. (Color online) Neutron ki- 
netic energy-differential cross section for the 
68 Ni(7*,n) 67 Ni channel. The experimental 
data (black points), the total neutron kinetic- 
energy fit (solid black line), the statistical de- 
cay component (dotted blue line) and the non- 
statistical decay component (dashed red line) 
are shown. 



performed on the reaction-probability level, the subtraction of the nuclear contribution requires 
the rescaling of the cross section measured on the carbon target for direct use on the lead data, 
using a semi-empirical scaling function based on the radii of the target nuclei. 

Before considering the curves obtained from the fitting procedure, one clearly sees that the 
experimental distribution extends significantly beyond neutron energies at which the 2n decay 
channel is open. Since for 68 Ni the In and 2n thresholds are located at 7.792 MeV and 13.60 MeV, 



respectively, the 2n channel opens at 5.81 MeV above the In threshold. It is therefore not 
surprising that the data above 5.81 MeV cannot be described by the statistical decay simulation 
alone, shown as the dotted blue line in Fig.[6j Since a different decay mode is required to describe 
the high-energy portion of the kinetic energies, the direct neutron decay will be investigated. As 
a rough assumption, only the direct decay into the ground state of the A-l nucleus is considered. 
The decay is easily simulated, since the neutron carries away the entire excitation energy minus 
its binding energy. The simulation of the detector response remains the same as in the case of the 
statistical decay. The obtained simulation data is shown in Fig. \6\ as the dashed red line, mainly 
describing the experimental data above 7 MeV. The sum of both contributions is represented 
by the solid black line, describing well the entire differential cross section distribution. A x 2 
minimization of the sum of the statistical and direct decay contributions leads to a direct decay 
branching ratio of 24(4)%, which is in good agreement with values obtained for other nuclei in 
this mass region [19j . 

5. El strength in 68 Ni 

The excitation energy-differential cross section of 68 Ni is reconstructed using Eq. [1] for the In 
and 2n data. Since the 3n threshold is located at 22.55 MeV, the analysis of the In and 2n 
channels covers most of the energy range of the GDR and of any low-lying strength. Using 
the procedure described in section [31 the main component of the trial input function consists 
of a Breit-Wigner distribution [1] used to describe the GDR. Since it is often more convenient 
to express the cross section in terms of El energy-weighted sum-rule strength, a fit parameter 
S ewsr is used to describe the fraction of the exhausted strength, in addition to the mean E m and 
width r parameters. In addition to the GDR function, five individual square bins, with widths 
ranging from 0.5 MeV at the In threshold up to 2.8 MeV at an excitation energy of 13.7 MeV, 
were added to describe possible deviations from a Breit-Wigner shaped El distribution in the 
In channel. The bin width was adjusted to take into account the resolution of the excitation 
energy, based on the statistical decay simulations described in section Before starting the fit 
procedure, the previously described strength distribution is converted to the proper Coulomb 
excitation differential cross section using the virtual photon spectra, as illustrated by Fig. 

The result of the fitting procedure is displayed in Fig. [7] for the measured Coulomb excitation 
cross section of 68 Ni. The fit of the GDR distribution was performed on the 2n data, since the 
eventual low-lying El strength is not expected to influence this neutron channel in any way. 
The x 2 minimization leads to E m = 18.1(5) MeV, T = 6.1(5) MeV and an El energy-weighted 
sum-rule strength Sewsr = 97(9)%, results that are in good agreement with systematic values, 
such as those proposed by Junghans et al. [20], predicting the GDR of 68 Ni at 18.17 MeV with 
a width of 5.17 MeV. Due to the width of the GDR and its proximity to the 2n threshold, a 
substantial amount of cross section belonging to the GDR is expected to reside in the In channel. 
The strong distortion induced by the experimental response however shifts this strength to lower 
energies, as shown by the blue dotted line in Fig. also clearly showing an excess of El strength 
below the GDR. It should also be noted that the fit procedure is carried out simultaneously on 
several observables, in order to obtain a more global description of the experimental data than 
by fitting each observable individually. In the present case, the fit is performed on the excitation 
energy, the neutron kinetic energy and their longitudinal momentum, both for the In and 2n 
decay channels, and additionally on the sum of the kinetic energies of the neutrons for the 2n 
channel. For instance, the fits presented in Fig. [6] and [7] were obtained simultaneously. 

A preliminary analysis of the additional El strength shows a Gaussian peak at 10.4(4) MeV 
with a width of a = 0.8(3) MeV. These parameters are in excellent agreement with the PDR 
measured by (7*, 7') virtual photon scattering [6], reported at 11 MeV with a width of less than 
1 MeV. In the present case, the peak integral corresponds to 4.1(1.9)% of the El energy-weighted 
sum-rule strength. This is again in good agreement with the value of Wieland et al. of 5% [6] 
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Figure 7. (Color online) Preliminary Coulomb excitation cross section of 68 Ni. The measured 
In (solid black points) and 2n (open red circles) data is shown, as well as the fit result for the In 
(solid black line) and 2n (dashed red line) channels. The dotted blue line shows the contribution 
of the GDR in the In data. 

obtained with a gamma decay branching ratio of 4% for the entire PDR region. 
6. Summary and outlook 

The analysis of the El strength in 68 Ni has been described, along with two important steps 
required to quantitatively extract its low- lying part. The first step consisted in evaluating the 
detector response of the R 3 B-LAND setup, resulting in the calculation of response matrices, 
allowing comparisons between theoretical strength distributions with the measured data. The 
second step was the measurement of the branching ratio for the non-statistical decay of the 
compound nucleus, also improving the quality of the fit of convoluted input distributions to the 
data. Based on this information, a fit of a trial El strength function, consisting of a Breit- 
Wigner function and of five independent bins spanning the In decay channel, was carried out 
and revealed a GDR in good agreement with systematics, as well as a non-negligible amount of 
low-lying strength. This preliminary additional cross section showed a peak structure exhausting 
4.1(1.9)% of the El energy-weighted sum rule. 

The observed PDR strength will provide a further constraint of the symmetry-energy 
parameters of the equation-of-state, similarly to Carbone et al. for the virtual photon scattering 
data on 68 Ni. This task will be performed using not only the El strength in the PDR region, 
but also using the dipole polar izability measured in the energy range starting at the neutron 
threshold and ending beyond the GDR. Since not only 68 Ni was measured in this experiment, 
but also more exotic species such as 70 Ni and 72 Ni, a similar analysis will be applied to these 
isotopes, allowing a verification of the behavior of the symmetry energy at different asymmetry 
values. 
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